Equation of state of the hot dense matter in a multi-phase transport model 
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Within the framework of a multi-phase transport model, we study the equation of state and 
pressure anisotropy of the hot dense matter produced in central relativistic heavy ion collisions. 
Both are found to depend on the hadronization scheme and scattering cross sections used in the 
model. Furthermore, only partial thermalization is achieved in the produced matter as a result of 
its fast expansion. 
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I. INTRODUCTION 

Many experimental results have been obtained from 
the collisions of heavy nuclei at the Relativistic Heavy 
Ion Collider (RHIC). To understand these results, both 
macroscopic hydrodynamic models [EJl, S 01 and mi- 
croscopic transport models [1, 0, 0, B Q have been used. 
Assuming that a thermalized quark-gluon plasma is pro- 
duced very early during the collisions, the ideal hydrody- 
namic model without viscosity is able to reproduce mea- 
sured transverse momentum spectra of various hadrons 
and their large elliptic flows. These experimental data 
can also be described by transport models that include a 
partonic phase in which partons undergo scatterings with 
cross sections much larger than those given by the pertur- 
bative QCD. Although the transport model approaches 
the ideal hydrodynamic model when the mean free path 
of a particle is much smaller than both the inter-particle 
distances and the size of the colliding system, it is not 
clear if these conditions are satisfied in relativistic heavy 
ion collisions. Studies have thus been carried out to un- 
derstand how results from microscopic transport mod- 
els differ from the macroscopic hydrodynamic models, 
particularly for the elliptic flow developed in relativistic 
heavy ion collisions [1, [l^l- Other questions which are 
also of great interest are the equation of state (EOS) of 
and the degree of thermalization in the hot dense matter 
described by these models. 

In the present paper, we study the EOS and pres- 
sure anisotropy of the hot dense matter produced in 
central relativistic heavy ion collisions in A Multi-Phase 
Transport (AMPT) model. The results indicate that 
both the EOS and pressure anisotropy depend on the 
hadronization scheme and scattering cross sections used 
in the model. Furthermore, only partial thermalization 
is achieved in the produced matter as a result of its fast 
expansion. 

This paper is organized as follows. In Sec. HIl the 
AMPT model is briefly reviewed. Results on the EOS 
in the AMPT model are shown in Sec. IIIII while the 
pressure anisotropy is studied in Sec. [iVl Finally, a brief 
summary and discussions are given in Sec. |Vl 



II. THE AMPT MODEL 

The AMPT model is a hybrid model [ll|, [H, [H, [11. 
It makes use of different models for different stages of 
relativistic heavy ion collisions. The initial conditions 
are taken from the Heavy Ion Jet Interaction Genera- 
tor (HIJING) [H]. In the default AMPT model, the 
partonic matter consists of only mini-jet partons from 
HIJING. It is different in the string melting scenario 
in which the hadrons generated by HIJING are dissoci- 
ated according to their valence quark structures and the 
resulting partonic matter is thus much denser. The evo- 
lution of the partonic matter is simulated with Zhang's 
Parton Cascade (ZPC) [H]- At the moment, only elastic 
parton scatterings are included in ZPC with the cross 
sections regulated by a screening mass. After partons 
stop scattering, they are converted into hadrons. In the 
default AMPT model, partons are combined first with 
their parent strings which then fragment according to 
the Lund model as implemented in PYTHIA [l^. In 
the AMPT model with string melting, closest partons 
are converted into hadrons via the coalescence model 
[13, [3 ■ The resulting hadronic system evolves accord- 
ing to A Relativistic Transport (ART) model [H, [l^]. 
The default AMPT model can give a good description 
of particle distributions [ll,[ll| while the AMPT model 
with string melting is needed to describe the elliptic flow 
[2l| and HBT radii [l^. Other observables have also 
been studied within the framework of the AMPT model, 
and these include phi meson production [23| , higher-order 
flows [13], pseudo-rapidity and svstem-si ze |26l. [27| 
dependence of flows, charmonium production |28l . |29| as 
well as charm [s^] and strange [3l| hadron flows. 



III. EQUATION OF STATE IN THE AMPT 
MODEL 



The equation of state of a matter can be described 
by the relation between its pressure and energy density. 
Since only particles contribute to the pressure and energy 
density in the AMPT model, the pressure and energy 
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density can be extracted from the energy-momentum ten- 
sor 



(1) 



where f{r,p) is the particle phase-space distribution 
function. When the latter is sampled by point parti- 
cles, the energy-momentum tensor can be calculated by 
averaging over particles and events in a volume i.e., 
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In Cartesian coordinates where the 2;-axis points along 
the colliding beam direction, the pressure components 
are related to the energy-momentum tensor by = T^^, 
Py = T^^, P^ = T"^^, and the energy density is given by 
£ = roo. In central collisions, because of the cylindrical 



P, 



the trans- 
P 



VI- 
For a 



symmetry around the beam axis, Px = 
verse pressure can be defined to be Pt = \{Px 
The longitudinal pressure Ph is the same as Pz- 
system in thermal equilibrium, its pressure is isotropic, 
i.e., Pt = Pl = P- Otherwise, we define the total pres- 
sure as P = \{Px + Py + Pz)- The equation of state of 
the system is then described by P/e as a function of e. In 
the following, calculations will be carried out for central 
Au-|-Au collisions at ^/Tnn = 200 GeV. In particular, 
we consider the central cell of produced matter, specified 
by the space-time rapidity 77^ and the transverse radial 
distance r. Only active particles, i.e., particles that have 
been formed but not yet frozen out, will be counted and 
contributions from free streaming particles are not in- 
cluded. 
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FIG. 1: (Color online) Equation of state in the AMPT model. 

Fig. [T] shows the equation of state extracted from 
the default AMPT model as well as the AMPT with 



string melting for three different parton scattering cross 
sections. In the initial partonic phase, the default 
AMPT gives a slightly higher P/e than the string melt- 
ing AMPT because the partonic phase in the default 
model mainly consists of massless gluons while that in the 
string melting model is dominated by light quarks with 
non-zero quark masses. Also, there are fewer partons in 
the default model than in the string melting model. Fur- 
thermore, hadronization happens earlier in the default 
model when e is around a few GeV while in the string 
melting model hadronization starts when e is an order of 
magnitude smaller. The partonic scattering cross section 
also affects the equation of state and the transition be- 
tween the partonic phase and the hadronic phase. When 
the cross section in the string melting model increases 
from 3 mb, to 6 mb, and then to 10 mb, the equation of 
state becomes harder and reaches into regions of lower en- 
ergy density. In the hadronic stage, the default model has 
a harder equation of state compared to that in the string 
melting model. The two models thus give very differ- 
ent pressures for the same energy density in the hadronic 
stage. 
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FIG. 2: (Color online) Equation of state in the partonic phase 
and the hadronic phase as well as combined partonic and 
hadronic matters from the default AMPT model. Arrows 
indicate directions of increasing proper time. The circles on 
the solid curve are data points taken at t = 1, 2, 3, . . . fm/c. 



The equation of state in the default model shows a 
sharp transition as a result of combining the equation 
of state of the partonic phase with that of the hadronic 
phase as shown in Fig. O In the default model, hadron 
production is related to the parent strings. The forma- 
tion proper time is the average parton freeze-out time 
plus 0.7 fm/c. As a result, hadron production is ear- 
lier for strings without final partonic interactions com- 
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pared with those that involve final partonic interactions. 
These early hadrons lead to the increase in hadron en- 
ergy density from about 3 x 10^^ GeV/fni'' at 0.8 fm/c 
when the system is still in the parton phase. Hadrons 
from strings with final state interactions lead to a second 
increase from about 2 GeV/inT' at 1.8 fm/c as seen in 
Fig. [21 The sudden change in the equation of state is 
due to the transformation of the string energy, which is 
not included in the calculations of the energy density and 
pressure, into hadrons during string fragmentation. This 
process completes at about 3 fm/c when the hadronic 
stage takes over. 



nate the energy density at about 13 fm/c. It is seen from 
Fig.[2]and Fig.[3]that the equation of state of the hadronic 
phase may have two different values for the pressure at 
a fixed energy density. This feature reflects the fact that 
the hadronic phase experiences the formation stage and 
expansion stage in the collisions, and the two stages have 
different particle compositions and/or momentum distri- 
butions and hence different pressures. 



IV. PRESSURE ANISOTROPY IN THE AMPT 
MODEL 
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FIG. 3: (Color online) Same as Fig.[2]from the AMPT model 
with string melting. 

The equation of state of the partonic phase as well as 
that of the hadronic phase in the string melting model 
with a parton scattering cross section of 10 mb are shown 
in Fig. [3] Unlike the default model in which gluons dom- 
inate the partonic stage, only quark degrees of freedom 
are included in the string melting model. As light quarks 
move out of the central cell, the relative abundance of 
strange quarks in the central cell increases, leading to a 
decrease of P/e as e decreases. The system stays in the 
partonic stage until about 7 fm/c when hadrons begin 
to affect the equation of state. The hadron energy den- 
sity at this time is about 10~^ GeV/fm^. As hadrons are 
formed from coalescing partons, their contribution to the 
energy increases. Although the energy density decreases 
as a result of the expansion of the system, the pressure 
drops even faster during hadronization. The equation of 
state of the whole system thus gradually changes from 
the partonic equation of state to the hadronic equation 
of state. Hadrons begin to play an important role when 
the hadron energy density reaches its maximum value of 
about 10^^ GeV/fm^ at about 10 fm/c and they domi- 
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FIG. 4: (Color online) Proper time evolution of the pressure 
anisotropy parameter in the AMPT model. 

In the above calculations, the energy density and pres- 
sure are obtained from the energy- momentum tensor, and 
they arc not necessarily for a matter in thermal equilib- 
rium. In order to check if the matter formed in relativis- 
tic heavy ion collisions is thermalizcd, we show in Fig. [4] 
the proper time evolution of the pressure anisotropy pa- 
rameter Pl/Pt defined by the ratio of the longitudinal 
to the transverse pressures. In the default model both 
the number of partons and the number of collisions are 
small during the early stage, and the anisotropy param- 
eter starts to increase quickly only when hadrons come 
into play at about 2 fm/c. In contrast, in the string melt- 
ing model there is a quick increase in Pl/Pt and by 1 
fm/c its value is already over 0.5. With increasing parton 
cross section, the system becomes more equilibrated and 
Pl/Pt becomes also larger. It is interesting to note that 
Pl/Pt crosses 1 at some time. However, it does not stay 
at 1 for any significant period of time. In the string melt- 
ing model, this happens after about 6 fm/c where there 
is a change from the gradual saturation of Pl/Pt to a 
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rapid increase, while in the default model, the change 
happens at about 10 fm/c. The proper time evolution 
of the energy density also shows a change at about the 
same time when the effect of transverse expansion sets in. 
After this time, particles with large transverse momen- 
tum move out of the central cell, leading to an increase 
in Pl/Pt- Therefore, to have Pl/Pt = 1 temporarily 
docs not indicate full thcrmalization during the evolu- 
tion of the system. Furthermore, Pl/Pt in the central 
cell can even go beyond 1. If only the initial stage is con- 
sidered, the anisotropy parameter can only reach a satu- 
ration value of about 0.8, implying that only partial thcr- 
malization is achieved in the central cell. We note that 
the present results are consistent with those obtained in 
Ref. with large two-body parton cross sections where 
the momentum anisotropy has been shown for the initial 
stage of central Au+Au collisions. It is also interesting to 
note that while larger cross section makes initial pressure 
anisotropy closer to 1, it also leads to delays both in the 
transverse expansion and in the rise of resulting pressure 
anisotropy. This is what underlies the crossing of the 6 
mb and 10 mb curves in Fig. 21 
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FIG. 5: (Color online) Time evolution of pressure anisotropy 
in the string melting model for two different space-time ra- 
pidity bin widths in the collider frame. 

For an accurate evaluation of the pressure anisotropy, 
momentum in the local rest frame needs to be used even 
for a space-time rapidity bin within |77s| < 0.5. Specif- 
ically, the local momentum p'^ can be calculated from 
energy E and momentum by 



with the proper time t = \Jt^ — [3^ . [ssj . The evolu- 
tion of the anisotropy parameter evaluated with the en- 
ergy and momentum in the lab or collider frame is shown 



in Fig. [5] The results depend strongly on the space-time 
rapidity bin width. However, if the local energy and mo- 
mentum are used, i.e., in the comoving frame, the de- 
pendence on the space-time rapidity width disappears as 
shown in Fig. [6l This indicates that the longitudinal flow 
is boost invariant around the mid-spacc-time rapidity. 
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FIG. 6: (Color online) Same as Fig. [5] but for proper time 
evolution in the comoving frame. 



V. SUMMARY AND DISCUSSIONS 

We have studied the equation of state in the AMPT 
model and found it to depend on both the hadronization 
scheme and the parton scattering cross sections. In the 
default model, the equation of state is much softer than 
that in the string melting model over a wide range of 
energy density as a result of early hadronization. In the 
latter case, hardness of the equation of state increases 
with increasing parton cross section. In both scenar- 
ios, the produced hot dense matter is not in full thermal 
equilibrium as the pressure isotropy deviates apprecia- 
bly from unity during most of its evolution as a result of 
fast expansion. The non-equilibrium nature of the equa- 
tion of state in the AMPT model thus differs from that 
commonly used in hydrodynamic models and also that 
obtained from lattice Quantum-Chromodynamics simu- 
lations [131 ■ There are recent attempts of improving the 
current description in the AMPT model by imposing 
hadronization at fixed time [3^1 ■ This certainly can make 
the equation of state at later time softer. Mean field can 
also modify the equation of state [s^. Also, an impor- 
tant aspect of the AMPT model is that the partonic 
phase includes only elastic scatterings. In a recent par- 
ton cascade model that includes parton number changing 
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processes @, the partonic matter is able to reach pres- 
sure isotropy. Parton inelastic scatterings are therefore 
important in bringing the system closer to ideal hydrody- 
namic solutions, which require at least pressure isotropy 
to be reached and maintained (STj . Furthermore, plasma 
instabilities may help the isotropization process as well 
[ssl . [39I . [40I . |41| |. The consequences of including these 
mechanisms in the AMPT model on both the equation 
of state of produced hot dense matter and other physical 
observables deserve further investigations. 
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